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I.  INTRODUCTION 


In  this  report  we  present  some  results  of  our  investigations  of 
structural  phase  transformations  in  media  forced  to  experience  time 
dependent  density  changes. 

Our  ultimate  goals  are  in  predicting  strain  and  temperature  ampli- 
tudes in  army  relevant  materials  which  undergo  structural  phase  trans- 
formations. While  those  ultimate  goals  have  not  yet  been  achievea, 
preliminary  results  give  promise  of  some  interesting  physical  interpre- 
tations, and  novel  military  applications.  We  shall  report  on  those 
preliminary  results,  and  their  possible  applications,  here. 

Phase  transformation  physics  has  progressed  from  a mature  and,  to 
some,  stale  branch  of  thermodynamics,  to  the  very  lively  study  area 
of  propagating  correlated  states1’2.  Those  propagating  correlated 
states  can  be  described  by  constant  shape  solutions  (called  solitons3) 
of  a nonlinear  equation  which  models  the  interatomic  potential  and 
statistical  mechanics  of  the  transforming  medium. 

The  approach  taken  here  will  be  that  of  the  middle  ground.  Rather 
than  follow  the  soliton  approach,  which  has  not  yet  been  applied  to 
media  of  dynamically  varying  geometry  (i.e.  time  dependent  strain), 
we  will  proceed  by  incorporating  some  of  the  more  modern  aspects  of 
phase  change  physics  into  a continuum  mechanics  framework. 

Our  interest  in  phase  change  phenomena  stems  from  optical  fuze  train 
considerations.  Picatinny  Arsenal  has  an  interest  in  developing  an 
optically  excited  detonator  containing  only  secondary  explosives  (for 
safety  reasons).  A typical  design  4 of  an  optically  excited  detonator 
is  shown  in  Fig  1. 


In  that  design  the  thin  metallic  layer  (e.g.  aluminum)  absorbs  the 
burst  (approximately  25  nanosecond  pulse  width)  of  optical  energy, 
gets  heated  beyond  its  vaporization  temperature,  and  expands.  That 
expansion  causes  a compressive  shock  wave  to  enter  the  explosive 
resulting  in  initiation  of  detonation. 

Unfortunately,  the  power  requirements  of  the  optical  detonator 
system,  illustrated  in  Fiq.  i,  for  an  explosive  such  as  RbX,  are 
approximately  4,5  4J/cm2  (25  nanosecond  pulse  width),  making  the  design 
of  a small  (and  thus  small  system  mobile)  laser  initiation  system 
difficult. 


Consider  the  thin  metallic  layer  shown  in  Fig.  1 for  the  case  of 
optical  energy  deposition  small  compared  to  the  metal's  vaporization 
energy.  In  this  case, ..the  pressure  induced  into  the  explosive  is 
proportional  to  the  Gruneisen  parameter,  r,  of  the  metal  ( r £ 9p/3e 
with  e the  thermal  energy  density).  Imagine  the  layer  to  be  10~4  cm 
thick,  1 cm  in  diameter,  and  to  possess  a typical  r of  2.  For  an 
absorbed  fluence  of  10"  2 J/cm2  (105  ergs/cm2)  the  cgs  pressure  generated 
in  the  metal  is  given  by* 


p = re  = 2 — cgs  = 2.6  x 10*  , (1) 

J (l)2  1C’4  cm 

or  P = 2.6  kbar  (1  bar  =106  dynes/cm2).  For  a specific  heat  of  0.215 
cal/gm°K  (0.90  J/gm°K),  corresponding  to  aluminum,  the  absorbed  fluence 
would  result  in  a temperature  rise  of 


AT 


10-2 

0.90 


2.7  J (l)2  10'4 


- 1 


°K  « 52°K, 


(2) 


*The  value  of  2.6  kbar  could  never  be  achieved  in  practice  under  the 
stated  conditions  and  25  x 10"9  sec-  optical  pulse  width.  This  is 
because  the  approximate  acoustic  transit  time  through  the  layer  (10~4 
cm/5  x 105  cm/sec)  is  10~2  times  that  of  the  laser  pulse  width  so  that 
the  generated  pressure  would  be  "leaked"  away  acoustically  before  the 
energy  deposition  process  is  completed  ( in  practice  one  might  expect, 
with  other  considerations  neglected,  a pressure  of  10“  2 x 2.6  kbar  or 
26  bar).  The  example  is  given  in  order  to  numerically  lead  into  phase 
transition  considerations,  rather  than  to  understand  the  detailed  physics 
of  laser  generated  pressure. 
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where  we  have  used  a density  of  2.7  gm/cc  (aluminum)in  the  calculation. 

The  2.G  kbar  pressure  of  Eq.  (1),  when  associated  with  a pressure 
pulse  of  less  than  10'6  sec  duration,  is  incapable  6 of  initiating 
detonation  in  secondary  explosives.  The  minimum  pressure  required  is 
of  the  order  of  20  kbar. 


' 


In  Fig.  2 we  show  Cruneisen  parameter  for  TiNi  as  a function  of 
temperature.  The  curve  was  obtained  7 by  measuring  the  thermal  expansion 

coefficient  while  the  TiNi  alloy  was 
undergoing  a structural  transformation 
at  approximately  60°C  and  atmospheric 
pressure  (during  a cooling  cycle). 
Clearly  a Griineisen  parameter  value 
associated  with  the  transition  peak 
in  Fig.  2,  when  substituted  into 
Eq.  (1),  would  meet  our  minimum 
pressure  requirements  for  initiation 
of  detonation  in  secondary  explosives. 
The  difficulty  is  that  the  large  r 
values  occur  over  a temperature  range 
of  approximately  5°C  while,  from  Eq.  (2) 
the  energy  deposition  is  expected  to 
result  in  a temperature  change*  of 
approximately  bO°K.  The  implication 
being  that  the  large  r values  are 
only  effective  over  a small  fraction 
of  the  energy  deposition  process. 
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Figure  2.  Griineisen  parameter 
vs  temperature  for  TiNi. 

Pace  and  Saunders7. 


From 


The  discussion  contained  in  the  above  paragraph  illustrates  the 
essentials  of  our  first  involvement  with  structural  phase  transformation 
physics  for  shock  wave  related  fuze  train  9 applications.  Clearly, 
while  there  is  promise  for  an  enhancement  effect  for  small  energy 
depositions  (so  that  the  resulting  temperature  change  does  not  carry 
one  beyond  the  transition  region),  the  effects  discussed  so  far  do 
not  offerany  promise  of  say  enhancing  2 kbar  to  20  kbar. 

There  is  another  approach  to  be  taken.  That  approach  is,  however, 
much  more  difficult,  and  requires  a deep  understanding  of  phase  trans- 
formation physics.  Materials  such  as  TiNi,  when  they  undergo  a phase 


* The  specific  heat  also  undergoes  large  variations  in  magnitude  over  a 
small  temperature  range  in  the  region  of  the  transition  temperature8. 
Because  those  variations  occur  over  so  narrow  a temperature  range,  the 
discussion  based  upon  the  50°K  result  of  Eq.  (2)  remains  valid. 
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transformation  either  give  up  energy  to,  or  absorb  energy  from  the 
surroundings.  That  energy  is  known  as  the  latent  heat  of  transformation, 
and  for  T i N i (51  atomic  percent  Ni ) the  latent  heat  is  5.8  Cal /gin 
(exothermic  on  cooling)11’10.  Now  the  question  becomes  can  we  pump 
(i.e.  force)  the  transformation  with  a shock  wave  in  such  a way  that 
the  liberated  heat  of  transformation  serves  to  reinforce  the  shock 
that  served  to  trigger  the  transformation.  A similar  situation  occurs 
with  the  detonation  process  in  explosives  - except  in  that  case  the 
integrated  heat  of  transformation  is  of  the  order  of  103Cal/gm. 

Figure  3 illustrates  what  we  have  in  mind.  Energy  density  and 
pressure  have  the  same  units.  Thus  5.8  Cal/gm,  for  a medium  of  6.5 
gm/cc  (NiTi),  is  equivalent  to  a thermal  pressure  of 

5.3  ~ = 38  = 160  ~ = 1.6  kbar.  (3) 

gm  cm3  cm3 

Consider  a slab  of  NiTi  alloy,  and  suppose  that  a pressure  pulse  of 
amplitude  0.5  £ 1.6  kbar  is  traveling  to  the  right.  In  the  differential 
volume  A assume  that  the  pressure  pulse  triggers  the  exothermic  trans- 
formation. The  result  of  that  triggered  transformation  could  be  a 
transformation  induced  pressure  pulse  of  0.5®  1-6  kbar  traveling  to 
the  left,  and  a similar  pulse  traveling  to  the  right.  If  the  correct 
conditions  prevail  the  two  right  propagating  pulses  (the  original 
pulse  plus  the  transformation  pulse)  could  reinforce,  resulting  in  a 
net  pressure  pulse  of  1.6  kbar.  An  adjacent  differential  volume  could 
give  similar  reinforcement  - the  net  result  being  a pressure  pulse  of 
avalanching  amplitude. 

There  are  a number  of  conditions  which  must  be  fulfilled  in  order 
for  the  avalanching  process  to  proceed.  While  we  will  not  now  dwell 
on  those  conditions,  let  us  simply  mention  that  transformation  induction 
time  and  the  post  transformation  thermal  equilibrium  relaxation  time 
are  physical  parameters  which  are  important  to  the  avalanching  conditions. 

There  are  a number  of  other  potential  military  applications  of 
structural  phase  transformations.  One  such  application  involves  the 
130  kbar  transition  n»12  in  iron  alloys  of  low  silicon  and  iron  content. 
That  transformation  offers  the  possibility  of  tailoring  armor  (or 
impacting  fragments)  so  as  to  control  free  surface  spallation  for 
offensive  or  defensive  purposes;  by  controlling  a transformation 
relaxation  parameter  one  can  apparently  control11  the  sharpness  of 
the  unloading  (rarefaction)  shock  and  consequently,  to  some  degree, 
control  spall  behavior.  We  shall  comment  further  upon  such  applications, 
and  the  relevant  transformation  physics,  within  the  main  body  of 
this  report. 
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II.  CONTINUUM  CONSIDERATIONS 

In  this  section  we  will  consider  some  continuum  aspects  of  static 
and  dynamic  phase  transitions.  We  will  begin  by  studying  a relatively 
poor,  but  instructive,  model  for  a transition  forced  by  a propagating 
strain.  Our  main  interest  will  be  in  further  evaluation  of  the 
possibility  of  reinforcement  avalanching  as  illustrated  in  Tig.  3. 

Consider  a one-dimensional  problem  characterized  by 

p0  - c(x,t)  0 = f ( x , t ) , (4) 

C ( X , t ) = cQ  + a(x,t)  [c,  - CQ  ] , (5) 

where  p is  mass  density,  u is  particle  displacement,  c is  an  elastic 
constant,  f dentotes  a viscosity  functional,  and  the  subscripts  zero 
and  one  denote  pre  and  post  transformation  states  respectively,  a can 
take  on  numerical  values  between  one  and  zero  and  is  meant  to  model 
the  physics  of  the  transformation.  In  writing  Eq.(4)  we  have  assumed 
an  ideal*  second  order  13  phase  transition  so  that  there  is  no  volume 
discontinutity  associated  with  the  transition,  and  with  the  consequence 
that  the  mass  density  reference  state  is  pQ  regardless  of  the  phase. 

Figure  4 gives  an  idea  of  the  intent  of  Eqs.  (4)  and  (b).  It  is 
known  from  experiment  14  that  tension  raises  the  transition  temperature 
(and  consequently  compression  lowers  the  transition  temperature)  in 
TiNi. 


*Ideal  second  order  transformations  are  characterized  by  volume  and 
entropy  changing  continuously  through  the  transition.  The  result  of 
that  continuity  is  zero  volume  change  and  zero  latent  heat  in  going 
through  the  transition.  Although  we  have  seen  that  the  NiTi  trans- 
formation does  possess  a non-zero  latent  heat  we  will  assume,  for 
the  purposes  of  modeling  with  Eq.  (4),  that  the  volume  changes 
continuously. 
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shift  in  transition  temperature). 


In  our  physics  picture  we  allow,  for  example,  a compression  pulse 
to  enter  the  solid  which  is  initial  1 1 y at  a temperature  T0.  The 
temperature  T0  is  assumed  to  be  only  slightly  below  the  zero  stress 
transition  temperature  Tc.  We  further  assume  that  the  compression  is 
strong  enough  to  change  the  transition  temperature  to  a value  Tc,. 

where  Tc  < TQ.  The  change  in  transition  temperature  with  stress  is 

very  dramatic;  for  bb  percent  (weight)  nickel  NiTi  al loy,  experiments14 
show  aT/aP  -v  lb°C/kbar.  Thus  a shock  amplitude  of  only  a few  kilo- 
bars  is  all  that  should  be  required  in  order  to  pump  a near  room  tem- 
perature transition. 

We  now  proceed  to  model  a in  the  form 


a(x,t) 


E (x , t ) t I 

V I . 


(6) 


where  e.  is  a critical  strain  energy  density,  t is  a relaxation  time 
characteristic  of  the  transition,  and  c is  the  impressed  strain  energy 
density.  Expanding  a to  first  order  in  (et)  gives 


3zu 

w 


i c 

i +j  £i  - 1 1 

et  32u 

!_  Uo  ) 

V *** 

= f(x,t; 


(7) 


Even  though  it  is  known  7,8  that  attenuation  (viscosity)  is  strong  in 
the  transition  region  we  set  f = 0;  we  have  not  yet  considered  possible 
physical  models  for  non  zero  f. 

2 

3 Li 

The  strain  energy  density,  e,  is  proportional  to—  , and 

i oX  j 

serves  to  make  Eq.  (7)  a formidable  partial  differential  equation. 

As  a consequence  we  are  forced  to  consider  two  different  approaches 
to  solutions.  The  first  is  analytic  and  very  approximate,  and  the 
second  is  via  coarse  finite  difference  techniques. 

A.  Approximation  I.  We  approximate  Eq.  (7)  with 


( 1 + et  ) 


0,  e = const. 


(8) 


Eq.  (8)  is,  neglecting  the  effect  of  finite  f,  a reasonable  approximation 
to  the  physics  for  the  following  two  conditions: 
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a.  The  frequency,  w,  associated  with  the  strain  is  so  large  that 
ut  >>1  holds.  Under  this  condition  we  can  argue  that  the  transforming 
medium  sees  an  rms  average  strain  rather  than  individual  peaks  and  troughs. 

b.  £ > £c,  so  that  the  strain  effect  on  the  transition  has  reached 
saturation  and  is  of  little  further  consequence.  This  condition  re- 
quires that  t < t in  order  that  the  a approximation  used  in  Eq.  (7) 
remains  valid. 


Neither  of  the  above  two  conditions  will  be  satisfied  by  the  subsequent 
algebraic  manipulations,  but  as  the  result  offers  a hint  of  relevant 
material  behavior  we  shall  proceed. 

Substitute  u(x,t)  = T(t)  x(*)  in  Eq.  (8)  with  the  result 


_ 1 

Ci+et)T 


-mV 


(9) 


where  -mV  is  the  separation  constant,  m is  a dimensionless  number 

o 

and  oj  is  a constant  with  units  of  angular  velocity.  T is  a dimension- 


less functional  while  x bas  the  dimensions  of  length. 
xm(x)  = )<„,  exp  -±  im%x 


Thus 


o J 


(10) 


and  the  T part  of  Eq. 


mV 

o 

' e2 


(9)  can  be  rewritten  as 
y T =0 


(11) 


with  y = (1+et). 

Equation  (11)  is  Airy's  equation  and  its  solutions  can  be  written 
in  the  form15 


Tm(t)  = A|nU1(Y,l)  + BmU2(Y,2), 


where 


( 1+  Bt), 


(12) 

(13) 


with  Uj  and  U2  being  linearly  independent  and  tabulated15.  The 
solution  u(x,t)  is  constructed  by  summing  the  product  Tm(t)xm(x)  over  m. 
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Let  us  now  evaluate  6 for  the  special  exaggerated  case  of  Cj=2c0. 


allow  £ = £ • 
c 

1 = 2co  ^ > 

.-i  . 

(14) 

= Cc 

^ i 

(m2u)  2t2)1/3 
0 

( l + t/t). 

(16) 

From  the  left  hand  side  of  Eg.  (9)  we  see 


d2T 

m 

_dt? 


m2io  2 


(16) 


so  that  if  we  identify  as  the  angular  frequency  of  the  applied 

disturbance  (evaluated  at  t=o)  there  is  justification*for  considering 
only  the  m=l  contribution  to  u(x,t).  We  thus  proceed  with  the  subscript 
m dropped. 


* The  amplitudes  Am  and  Bm  are  determined  by  the  orthogonality  of  the 

set  of  functions  of  Eq.  (10).  Let 

u(x,t)  = T.  T (t)  xn(x). 
n n n 

Multiply  both  sides  by  exp 
over  x 


Po  1 

± i mw  x,—  : , and  integrate 

0 v c0  ! 


/j*J  0 

e * V c^  u(x,t)dx  = 2n  E Tn(t)xm  6(n-m). 

- QO 

Tm(o^xm  = Jt\  Je  u(x,o)  dx. 


The  above  result  implies  that,  for  an  ideal  disturbance  which  at  t=o 
is  a pure  sinusoid  (of  angular  velocity  u>p)  filling  all  space,  the 
solutions  are  characterized  by  Aj  and  Bj  being  the  only  amplitudes 
which  do  not  vanish  identically.  While  that  ideal  case  can  never  be 
achieved  in  practice  (in  a typical  experiment  the  t=o  disturbance 
resides  only  within  the  boundaries  of  a transducer),  the  ideal  case 
serves  as  the  foundation  for  concentrating  our  attention  only  on  the 
m=l  component  of  the  solution. 


From  Eqs.  (10),  (12),  (13),  and  (15) 


u(0,0)  = 

AU.  (u  t)2/3,  1 1 
1 0 

+ BuJ(u)  t)2/3,  1 

2 [ o 

3u(0,t) 

at 

A (V)2/3  Uf 

T 

(V)2/3.  i j * 

0 

+ B (V)2/3  U2 

'<V)2/3,  i . 

T 

functions  have  the  properties15 

Uj  (0,1)  = 1,  U2  (0,1) 

= o. 

Uf  (0,1)  = 0,  U2  (0,1) 

= 1. 

(17) 


(18) 


(19) 

(20) 


We  now  make  the  "hydrodynamic"  approximation  of  (w  t)  <<  1 (which 

is  contrary  to  the  condition  (a)  which  follows  Eq.  (3)  ) so  that 
from  Eqs.  (17)  and  (19) 

u (0,0)  = A,  (21) 

and  from  Eqs.  (18)  and  (20) 


3u(0,t)  I . B (u  t)2/3.  (22) 

3t  | t 0 


It  has  already  been  assumed,  see  Eq.  (16),  that  the  disturbance 
may  be  characterized  at  (x,t)  = (0,0)  by 


au(0,t) 

at 


O 


'v.  a)  u(0,0). 

'V/  0 


(23) 


Eq.  (23)  is  a rise  time  statement  rather  than  a requirement  that  the 
applied  disturbance  be  a plane  wave  sinusoid.  Thus 


u (0,t)  = u (0,0) 


MY.l)  + (V)1/3  U2 ( Y,  1 ) , 


(24) 
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which  predicts  that  in  the  hydrodynamic  approximation  u(0,t)  has  the 
time  dependence  of  Uj(Y,l).  Figure  5 displays  that  time  behavior. 
From  Eq.  (15)  we  see  that,  for  t>0,  our  interest  resides  in  Y>0  for 


In  order  to  judge  the  importance  of  the  physics  modeled  by  Eq.  (8) 
we  must  compare  displacements  with  and  without  the  presence  of  the 
modeled  phase  transition.  That  comparison  can  be  achieved  by  assuming 


Figure  5.  U,(Y,1)  and  U2(Y,1)  vs  Y. 

2/3 

(a)  t)  ' corresponds  to  t=0  for  c,  = 2c  , 

0 2/3  0 

and  (2w0x)  corresponds  to  t=0  for 

co  = 2ci  ' 

U(0,t)=u(0,0)cosu>  t (25) 

o 

in  the  absence  of  a transition,  and  inquiring  as  to  the  value  of  t 
(which  we  call  t,)  at  which  Uj(Y,l)  achieves  it  first  node.  From 
Eq.  (15)  with  m=l,  and  Fig.  5, 

(V)2/3  (1+  li)  ■ 2 (26) 
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For 


L t)  = 0.1  we  thus  have 

o 

i + io  (^tj)  = |onr^3  ^ ti  = * (27) 

t compares  to  a time  of  tt/2w  = 1.57gj  for  the  node  to  be  achieved 

i 0 0 

in  the  absence  of  a transition.  That  difference  in  time  increases 

as  (w  t)  decreases. 

o 

The  above  analysis  shows  that  the  modeled  transition  serves  to 
accelerate  changes  in  the  displacement  u(o,t).  That  acceleration  may 
be  viewed  as  a rudimentary  form  of  avalanching.  For  the  cx  = 2cq  case 

that  avalanching  is  limited  by  the  oscillatory  behavior  of  U2(V,1) 
for  Y>0.  It  should  also  be  noted  that  the  result  of  Eq.  (27)  violates 
the  condition©  (which  follows  Eq.  (8)  ).  That  condition  was 
important  for  interpreting  the  physics  appropriate  to  the  model  of  Eq.  (8). 

One  can  also  consider  the  (exaggerated)  elastic  constant  cnange 
of  Cx  = cq/2.  Such  a change  can  be  obtained  by  requiring  (see  Fig.  4) 


T 

o 


rather  than 


ci  = 2co 


(28) 


(29) 


as  illustrated  in  the  two  examples  of  Fig.  4.  For  the  change 
c = cq/2  , Y and  u(o,t)  become16,  subject  to  the  same  assumptions  as 

in  the  discussion  of  ci  = 2cq, 


Y = (2u  x)2/3 
0 


(30) 


u(o,t)  = u(o,o) 


Ui(Y,l) 


(2u)ot)1/3  U2 ( Y , 1 ) 


(31) 


Equations  (30)  and  (31)  correspond  to  decreasing  Y for  increasing  t. 
Here  the  avalanching  effect  is  quite  dramatic  (see  Fig.  b)  in  that 
Ui(Y.l)  is  unbounded  for  Y<0.  As  with  Eq.  (24),  Eq.  (31)  suffers  In  that 
its  use  violates  conditions  © and<E)which  follow  Eq.  (8). 


13 


However,  regardless  of  the  violation  of  those  conditions,  the  above 
results  hint  at  a type  avalanching  behavior  consistent  with  Fig. (3). 


B.  Approximation  II.  We  now  will  attempt  to  apply  crude  finite  dif- 
ference techniques  to  Eq.  (7)  with  f = 0 . Our  objective  being  to  gain 
physical  insight  beyond  that  capable  with  Eq.  (8). 


In  Eq.  (7)  with  f =0  substitute  vQ2  i c0/Po,  v22  = Ci/p0, 


e s m , 


au  y 
'v  ax/ 


and  cc  = Mn2. 


M is  to  be  identified  as  an  elastic 


constant,  nc  a critical  strain  value,  V]  the  longitudinal  sound  velocity 
in  the  transformed  medium,  and  vQ  is  to  be  identified  as  the  pre- 
transformation longitudinal  sound  velocity.  The  result  of  the 
substitutions  is 


a2u  2 32u 

at7  " o ax7 


(32) 


We  now  employe  the  so-called17  characteristic  stretching  trans 
formation 

5 = x - Vt,  c 5 aVt.  (33) 


We  further  define 

»(5,C)»~.  (34) 

The  concept  associated  with  Eqs.  (33)  is  that  of  a disturbance 
characterized  by  an  almost  constant  phase  velocity  V;  that  is  we 
seek  solutions  of  Eq.  (32)  which  have  only  weak  c dependence  (and 
conversely  a strong  £ dependence).  This  is  equivalent  to  assuming 
a small  (and  neglecting  terms  in  a2).  If,  for  example,  the  non- 
linear term  in  Eq.  (32)  were  zero,  then  a solution  would  be  exp  (ike) 
with  k constant  and  a identically  zero.  In  other  words  we  seek 
solutions  which  are  only  slightly  removed  from  plane  waves. 


Substituting  Eqs.  (33)  and  (34)  into  Eq. 
neglected,  gives 


(32),  with  a2  terms 


2aV2f  + 
c 


V2)  + (V!2  - Vo2) 


(35) 
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where 


C,  h aVx . 

a 


(36) 


Straight  forward  finite  differencing  of  Eg.  (35)  results  in 


, <V  W ' * <«„•  Sn1 


J - 

Fl  - FlV!|V  r-m> 

^n+l.Sj  " 4'^n,SJ 

A 

. 

J 

where 


’ (37) 


r V2  - V-2  v 2_  V 2 

Fj  - 2_  , and  Jr 2 H 1 o 


2aV2 


2aV2canc2 


(38) 


6 and  a are  mesh  dimensions  and  are  defined  by 


?m+l  = <m  + 5 ’ and  Vl  = *n  + 


A. 


(39) 


Rewriting  Eq.  (37)  with  the  time  derivative  single  stepped  backwards 
gives 


y(n,m)  - y(n,m-l)  = [g 

t - G2cmy2(n,m)]  y (n+1 ,m)  + 

- [•: 

[ - G2tmy2(n,m)]  y(n,m). 

(40) 

where 

Gi  r -f  Fi  and 

u 

G2  -=  i F2  • 

(41) 

Thus,  solving  Eq.  (40)  for 

y(n+l,l)  - y(n,l) 

y(n+l,l)  - y(n,l)  = 

bi  - 

- y(n,o) 

(42) 

If  the  disturbance  enters  the  transforming  medium  from  the  left 
at  (x,t)  = (0,0),  then  the  t=0  strain,!^, vanishes  for  x>0.  Ue  thus 

9x 

Furthermore,  for  t = m 6 and  Cj  = a V6  Eq.  (42) 

d 


set  y(n,0)=0  for  n>0. 
becomes  for  n>0 
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’("♦'•'i  - r 1 - 

d c o 


-i 


(43) 


For  'V2  -o  n 2 and  m > 1,  assuming  vt2-v  2 > V2  -v  2 , Eq.  (43)  can 

< C d *0^0 

be  approximated  by 


f(n+l ,1 ) - 'F(n.l)  ^ 


Or 


T(n,l) 


, , '*,?(n,l)  v>  - Vo 

1 + m'T"?"  Z~r~ 


manc 


o 

2 - v 2 

V2" -v 


A \ 


f(n+l,l)  Y (n,l)  + 2a(| 


V2 


V2_  v 2 
0 


16/  manc2 


f(n , 1 ) + 

V2(Vl2  -v  2)1 

LF^vtj 


(44) 


(4b) 


The  second  term  on  the  right  of  Eq.  (45)  gives  a purely  coordinate 
system  dependent  linear  strain  contribution  to  the  space  dependent 
strain  growth.  The  third  term  on  the  right  contains  the  transformation 
physics(again  coordinate  system  dependent  because  of  the  presence  of  a) 
contribution  to  the  spatial  strain  growth.  The  third  term  is  cubic 
in  the  strain,  positive  for  C!  = 2cq,  and  negative  for  Cj  = c0/2. 

Equation  (45)  predicts,  for  t < t (which  allows  the  approximation 
of  Eq.  (44)  ),  strain  as  a function  of  position  for  a medium  undergoing  a 
strain  pumped  transformation.  It  is  consistent  with  spatial  avalanching 
in  that  it  contains  a cubic  strain  contribution  to  the  strain  growth. 

An  intelligent  choice  for  V could  be  an  experimentally  observed18 
velocity  for  propagation  of  a dynamic  transformation. 

The  two  continuum  approximations  considered  in  this  section  hint 
that  shock  amplification  can  occur  in  the  presence  of  a strain  pumped 
structural  transformation. 


III.  LATTICE  CONSIDERATIONS 

Solid-solid  phase  transformations  occur  via  two  possible  micro- 
scopic mechanisms.  Broadly  stated  the  mechanism  of  transfromation 
related  atomic  motion  may  be  classified  as  either  diffusive  or  diffusion- 
less. 

Diffusive  transformations, such  as  that  associated  with  the  formation 
of  pearl ite  in  low  carbon  content  iron  alloys,  are  slow  processes19. 
Conversely,  diffusionless  transformations,  such  as  that  associated 
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with  martensite  formation  in  quenched  steels19,  can  occur  quite 
rapidly  18>19.  The  near  room  temperature  transformation  in  almost 
stoichiometric  NiTi  is  of  the  diffusionless  variety10*8  and  is  often 
called  a "martensitic"  transition10*8. 


A rough  idea  of  the  maximum  transition  time  for  a transformation 
to  be  important  to  shock  wave  physics  can  be  obtained  by  considering 
the  planar  impact  between  two  half-centimeter  discs  of  identical 
material  as  shown  in  Fig.  6.  The  impact  initially  produces  a compres- 
sive shock  of  pressure  amplitude  p0VUs/2  (see  F’ig.  6 for  definition 
of  symbols),  and  that  pressure  is  reduced  to  zero  when20  the  re- 
flected shock,  of  pressure  amplitude  - p0VUs/2,  uncovers  a previously 
compressed  differential  volume.  A differential  volume  is  thus  at 
pressure  for  a time  t!  where  tj  <_  «. x / Us  Thus,  to  the  extent  that 
shock  pressure  triggers  the  transition,  the  transition  must  occur 
in  a time  less  than  tj.  For  0.5cm  and  Us  = 5xl05  cm/sec,  t,  <_ 
2ysec. 


Another  gross  criterion  for  shock  pumping  a transition,  with  con- 
sequent reinforcing  of  the  original  shock,  can  be  stated.  If  the 
transformation  can  propagate  at  a velocity  approximately  equal  to  or 
greater  than  the  shock  velocity,  then  shock  pumping  with  reinforce- 
ment is  a possibility.  If  the  shock  transformation  propagates  at  a 
velocity  considerably  less  than  the  shock  velocity,  the  transition' 
may  be  pumped,  but  reinforcement  probably  can  not  occur. 


It  is  of  interest  to  compare  diffusive  transformation  numbers 
with  the  above  velocity  criterion.  For  a concentration  n the  diffusion 
coefficient  D is  defined  by  Fick's  equation 


M = n A AA 

At  AX  [AX 


(46) 


For  a diffusive  displacement  mean  free  path  i and  corresponding  mean 
free  time  x,  Eq.  (46)  can  be  approximated  by 


(47) 


Thus,  for  a diffusion  velocity  VD  defined  by  VQ  = t/x,  we  have 


VD  - DM 

For  a solid  or  Mquid  the  relevant  mean  free  path  is  approximately 
a lattice  parameter  (3xl0_8cm),  while  for  a gas  minimum  t is  identi- 
fied as  the  mean  molecular  separation.  For  an  STP  gas  l % 3.5xl0*7cm. 
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Table  1 lists  the  three  different  ranges  of  the  diffusion  constants 
with  the  resulting  diffusion  velocities. 


Table  1 

Diffusion  Coefficient3  and  Velocity 

Substance 

D 

VD 

gas 

liquid 
sol  id 

- . 

0.2  cm2/sec 
10'5 
10" 1 2 

5. 7xl05  cm/ sec 

3xl02 

3xl0*5 

Diffusion  coefficients  obtained  from  American  Inst.  Phys.  Handbook, 
D.  E.  Gray,  editor  (McGraw-Hill,  New  York,  1972).  Third  edition. 

It  is  obvious,  from  Table  1,  that  shock  pumped  diffusive  trans- 
formations in  liquids  and  solids  are  not  expected  to  give  rise  to 
reinforcement*.  As  our  interest  is  in  solids  we  will  now  consider 
some  physics  of  diffusionless  transformations. 

Phonon  Formulation.  While  diffusive  processes  are  characterized  by 
random  atomic  (or  molecular)  displacement,  diffusionless  processes 
are  associated  with  coherent  (the  opposite  of  random)  displacements. 
The  displacements  are  coherent  because  there  is  a definite  phase 
relationship  between  the  displacive  transitions  of  neighboring  atoms. 

The  thermal  vibrations  of  the  atomic  constituents  of  a crystalline 
(lattice)  solid  can  be  described2 1 »22  by  a set  of  independent  harmonic 
oscillator  states.  The  degree  of  occupancy  ( and  the  corresponding 
vibration  amplitude)  of  each  state  is  a function  of  the  oscillator 


* Gases,  with  STP  acoustic  velocities  in  the  region  of  3.3xl04cm/sec, 
would  appear  to  be  prime  candidates  for  reinforcement.  Perhaps  that 
is  the  reason  why  explosive  detonation  in  gaseous  mixtures  is  so 
easy  to  achieve. 
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frequency  and  the  temperature.  The  smallest*  change  in  occupancy  of 
a given  state  is  called  a phonon  21*22.  The  elementary  (i.e.  differential) 
oscillations  of  each  phonon  in  a given  state  are  exactly  in  phase 
(i.e.  coherent),  and  summing  over  the  different  states  gives  the  net 
vibration  of  the  crystal  with  respect  to  its  center  of  mass. 

The  displacement  x of  a simple  harmonic  oscillator  of  mass  m 
obeys  the  differential  equation 

Ht*  + “2x  =0’  “2  = m ’ (49) 


where  u>=  2nf  with  f being  the  oscillator  frequency,  and  K is  the 
"spring  constant."  The  spring  constant  can  be  expressed  in  terms  of 
elastic  constants  which  are  in  turn  a function  of  crystal  structure  23,24 
(or  equivalently  a function  of  the  interatomic  potential),  u is  a 
function  of  the  state  variable,  and  for  an  ideal  solid  has  a functional 
dependence  as  shown  in  Fig.  7. 


★ 


A change  in  occupancy, 


, causes  a change  in  state  energy  via 


(AEl  = (AN)  hflq] 

k S ^ 

where  h is  Plank's  constant,  and  f (q)  is  the  oscillator  frequency 
for  the  vector  state  g (a  N)  = + 1 defines  the  addition  or 


subtraction  of  a 


single  phonon  with  respect  to  the  vector  state  of  q 

'V, 
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and  9 » 9„  for  some  critical 

vector  , where  Tc  > 0 • 


Figure  7.  Frequency  spectrum 
in  [lOOj  direction  for  an  ideal 
cubic  lattice  of  unit  cell  dimen- 
sion a. 


shock  wave  effects  in 
that  oj  > o as  T — ■*  Tc 

Tc,  and  a critical  state 

Such  behavior  is  related  to 
a structural  instability25 
(such  as  associated  with  a 
phase  transition)  and  is  not 
a property  of  an  ideal  lattice 
as  illustrated  in  Fig.  7.  A 
solid  having  the  above  property 
is  said  to  possess  a soft  mode25*26 
There  is  a point  of  view27  that 
all  phase  transformations , in- 
cluding the  liquid  < - ) gas 
transformation,  can  be  described 
in  terms  of  soft  mode  physics. 

The  "mode"  ufit)  is  said  to 
be  soft  because,  from  Eq.  (49)> 
when  oj — > 0 the  spring  constant  K 
also  goes  to  zero  (i.e.  the 
spring  becomes  soft).  As  the 
vibrational  restoring  force  is 
proportional  to  K,  K — » f)  implies 
that  the  displacement  does  not 
return  to  zero.  A displacement 
not  returning  to  zero  is  con- 
ceptually equivalent  to  a 
structural  transformation.  For 
subsequent  analysis  we  shall 


In  this  section  we  will  be  interested  in 
those  lattice  solids  which  have  the  property 

temperature 


employ  the  functional  form25 


-2(!)  - * (T  - Tc> 


+ B 


q 

'V, 


(50) 


which  exhibits  the  soft  mode  property. 

We  begin  by  considering  a modified  Fourier  transform  for  the 
macroscopic  displacement  u(x,t) 


a 

u(x,t)  = 1 f 

\T*  J 


o (q)  N (q,x,t)  uq  (t)  e iqx  dq,  (51) 
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where  G (q)  is  a geometry  dependent  (density  of  states)  factor, 
N(q,x,t)  is  the  previously  mentioned  occupancy  of  the  state  q,  and 
Uq(t)  is  the  displacement  amplitude.  In  principal  uq(t)  contains 

thermal  vibrations  (with  a soft  mode  contribution),  and  also  those 
displacements  which  constructively  sum  to  give  the  center  of  mass 
motion  necesary  for  successful  continuum  mechanics. 


Equation  (51)  describes  only  one  dimensional  physics.  Although 
structural  transformations  can  occur  on  a one-dimensional  lattice28, 
Eq.  (51)  is  limited  in  application  in  that  the  general  transformation 
is  explicitly  three  dimensional29*.  We  can,  however,  still  carry 
out  some  rather  interesting  and  meaningful  analysis. 


Operate  on  u(x,t)  of  Eq.  (51)  with  the  wave  operator  W(x,t), 


W(x,t) 


(52) 


uq(t)  + 


eiqx  dq.  (53) 


In  the  absence  of  interesting  effects  such  as  viscosity,  non- 
linearity, or  mode-softening,  we  would  expect  the  right  hand  side 
to  equal  zero  (corresponding  to  propagation  of  a simple  disturbance 
with  velocity  (co/p0)re  ).  Since,  however,  we  are  concerned  with  a 
particular  phonon  state  7^.  ( the  state  with  a softening  spring  constant) 

we  can  argue  that  to  lowest  order  the  physics  of  primary  interest 


* Structural  transformations  are  usually  discussed  in  terms  of  a 
discontinutity  in  the  three  dimensional  symmetry  transformation 
group  (the  set  of  transformations  which  leaves  the  structure, 
or  its  Hamiltonian  equivalent  unchanged). 
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resides  only  in  the  direction  qc/|qc|.  It  is  to  be  understood  that 
Eq.  (51),  and  the  consequent  analysis,  has  been  written  in  that  spirit. 

The  first  curly  bracket  on  the  right  in  Eq.  (53)  represents  second 
sound30(propagation  of  a distrubance  in  occupation  number).  Since 
pure  second  sound  is  not  accompanied  by  a mass  density  distrubance30, 
it  will  be  naively  assumed  that  the  contents  of  the  first  curly  bracket 
identically  sum  to  zero.  Any  first  curly  bracket  contribution  to 

Wu(x,t)  awaits  further  work. 


Now  consider  the  aN/ax  contribution  to  the  second  curly  bracket 
in  Eq.  (53). 


p / \ 9 N 
G(q)  a? 


Ug ( t ) e iqx  dq. 


(54) 


But  in  equilibrium 


N (q,x,t)  = 


exp 


w = w + Aw  , 

o 


+ AT, 


(55) 

(56) 


whereli  e h/2-n , k is  Boltzmann's  constant,  and  Aw  and  AT  are  shock 
related  changes  in  the  mode  angular  frequency  and  local  temperature. 
Using  Eqs.(56)  to  expand  Eq.  (55)  to  first  order  in  Aw  and  AT  gives 


+ 


x0e  xo 
ex°-l 


(57) 


where  xQ  E"fc  <WkTQ  and  Nq  = NQ(q)  is  the  pre-shock  occupation  number. 
Thus 


3N 

3X 


H 

3X 


w 3X 
o 


(58) 
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by 


In  phonon  language  the  Griineisen31  ,32  parameter  r (q)  is  defined 


(Aw(q)  ) = - r (q)uo  (q)  (An), 


(59) 


with  n be  the  local  strain  (n  = |~) . Substituting  Eqs.  (58)  and 
(59)  into  Eq.  (54)  gives 

|2  \1/2  co 


^*lu(x,t)^2  = 


Mn 


poTo 


il 

ax 


o 

/ 


q G (q)  x eX°  N 2 u (t)e'H*dq  + 
o o q 


iq  xr 


- 1 


p a/2  c 
- 1 -2- 

(3n| 

I11  1 Po 

|ax 

o q 


For  a mode  independent  Griineisen  parameter  (i.e.  r(q)  = r = constant) 
one  can  show33  that  (AT)/Tp  = - rn  with  the  result  that  ^Wu(x,t)^ 2=0. 
For  a transforming  solid  the  assumption  of  constant  r is  particularly 
poor;  in  the  introductory  section  (see  Fig.  2)  it  was  noted  that  r 
is  a strong  function  of  temperature  in  the  vicinity  of  the  transition 
temperature.  Current  thinking34,  at  least  for  ferroelectrics , is 
that  strong  temperature  dependence  is  related  to  the  behavior  of 
r(g)  near  £c.  If  we  thus  let 


r(q)  = r + (Ar)q  (61) 

in  Eq.  (60),  and  define 

<(Ar)q>i  J q G (q)  Ur)qxoeX°  Nq2  Uq(t)  e 1qx  dq,  (62) 

Eq.  (60)  reduces  to 


{Wu(x,t)>2 


^2  £T 

Po  3x 


(63) 


where  g2  plays  the  role  of  a thermoelastic  coupling  coefficient. 


g2  = - 


<(AF)q>. 

r 


(64) 
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The  terms  remaining 
into  what  we  shal 1 cal  1 


to  be  considered  in 
{Wu(x,t))1 


Eq- 


(53)  can  be  grouped 


{Wu(x,t))j 


G (q)  N (q,x,t)  e iqx© 


3N 

at  dt 


dq. 


(65) 


As  previously  mentioned,  each  discrete  q value  defines  an  independent 
harmonic  oscillator  state 


+ 1.  ^ + “>2(q.T)u  = 0,  (66) 

IF  * t dt  q 


where  allows  for  the  presence  of  damping25.  We  have  not  explicitly 

included  a strain  dependent  forcing  term  in  Eq.  (66).  The  idea  is  to 
implicitly  include  such  an  effect  via  the  strain  dependence  of  T and 
Tc  through  Eq.  (50).  Substituting  Eq.  (66)  into  Eq.  (65)  yields 


{Wu(x,t) } 


G (q)  N (q,x,t)  e 


iqx 


© 


.2  3N  1 

N at  ‘ 


xq  / dt 


} dq. 


(67) 


The  evaluation  of  Eq.  (67)  will  begin  with  a rather  drastic 

simplifying  assumption.  It  will  be  assumed  that  Eq.  (67)  vanishes 

except  near  q = q , and  that  near  q 
v C 


1_  „ 2 3N 

xq  ‘ N at  * 


(68) 
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Eq.  (68)  states  that  the  phonon  relaxation  effects  are  just  balanced 
by  the  dynamic  strain  induced  phonon  source  term  M.  To  claim  other 

than  Eq.  (68)  requires  more  detailed  physics  than  is  currently  avail- 
able to  us.  Thus,  upon  employing  Eq.  (50)  at  q = qc>  and  Eq.  (68), 

in  Eq.  (67)  gives 


(WU(x,t»i  ^ G(qc)  N(qc)  j q/  £ - A <T-TC)|  u, 


iq  X -iqx 
® { e + e 


(69) 


where  we  have  assumed*  that  G,N,  and  un  are  even  function  of  q near 

4c 

q = qc>  (Aq)  denotes  a bandwidth  centered  around  qc. 

The  analysis  which  follows  below  utilizes  {Wu(x,t)}. 

If  the  further  drastic  assumption** u = const.  = u is  made 

qc  ^ 

Eq.  (69)  reduces  to 


fc<c>j!r 

■ar[  cos  V* 

(70) 

2 \i/2 

<“»(!] 

G (qc)  N(qc)  A Uq  . 

(71) 

* The  assumption  is  trivial  for  G which  happens  to  be  a constant  for 
one  dimension  geometry  in  q space.  We  have  retained  G(q)  in  our 
calculations  as  a preliminary  to  eventual  multi-dimensional  calculations. 

**This  assumption  means  that  u represents  the  correlated  lattice 

qc 

displacement  which  results  in  a structural  transformation. 
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A.  Thermoelasticity  with  a Soft  Mode.  In  a previous  report35  it 
was  shown  that  thermoelasticity  in  the  presence  of  thermal  conductivity 
and  longitudinal  to  transverse  temperature  (phonon  occupation  number) 
relaxation  satisfies 


(1  + 


It)  <1+Ti  It)  T ‘ T + 


at' 


T„tlT 

0 0 1 

2 g p 
yo  o 


T3  It) 


V0TlT2 


(1  + T1  It} 


a2t 

lx7 


(72) 


where  v(x,t)  = , n(x,t)  = - 9u^-,t:)  , (73) 

(notice  the  change  in  sign  convention  for  n from  that  used  earlier 
in  this  report). 


aTy  . i_ 

at  = 


(74) 


with  Ty  and  Tz  being  the  temperatures  transverse  to  the  x direction 
of  strain  propagation  (T  =T),  and 


aT  , T h 

— * = 1 (j  . t ) + _2_o 
at  T,  Uy  pogo 


la  i pTx 

at  at 


*TC 


(75) 


which  g is  given  by  c /3  (with  c„  being  the  specific  heat  at  constant 
o /aT  \ v sT 

volume)  and  x)  is  the  thermal  conductivity  contribution35  to  x . 

4c 


at 


!L> 

at 


TC 


V*2  l>2\ 

2 ,ax2 


(76) 


In  Eqs.  (72)  and  (76)  vQ  = (Co/Po)  1/2 


27 


Before  proceeding  with  the  analysis  it  remains  to  find  the  re- 
lationship between  the  thermoelastic  coefficient  hQ  (of  Eqs.  (72) 

and  (75)  and  of  Eq.  (63).  A linear  thermoelastic  medium  has  a 

coefficient  hQ  defined  by 


(77) 


where  a is  stress.  Thus  the  strain,  corresponding  to  stress 
free  thermal  expansion  has  the  property 


3n 

3x" 


TE 


II 

3X 


(78) 


But,  from  Eq.  (63),  considering  time  independent  thermal  expansion 

with  n = - ISi 

3X 


- ~0  3ZU 


p 3X^ 
o 


_o  _3_n  _ 

p 3x 
0 


92 


3T 


Po  3X 


(79) 


with  the  result  that  h0  = 92  • Thus  combining  Eqs.  (63)  and  (70)  gives 
h 


(“#)  = - r - F<o 

0 


9l 

at  " at 


(80) 


Proceeding  now  with  the  analysis  by  operating  on  Eq.  (80)  with  j_. 

3X 


3 V 


0 3 1 32T 

) 

0 


» - P-*  ft  fxrr  f«u 


32T  ,nc 
3xat  ' 3x3t 


cos(qox)  + 


* <c  F(<»C) 


il  I[c 
3t  " 3t 


sin(qcx). 


(81) 
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OOlCvJ 


Now  operating  with  W of  Eq.  (52)  on  Eq.  (72),  and  substituting  from 
Eq.  (81),  gives 


' "T ' ^v7 (lt,>  5t>  Tt  (0i  * 

T h T a j\2T  _ + 

•F(qc)  44“  C0S  (qcx)  (1  + T1  ft}  nix  Iu7 
oyo  L 

i 

+ qc  F(qc)  o o i sin  (v°  (1+  T,  it  ' it  ' w " 

L J 


V 2 T.T,  3 a2T 

= -JLVJL  W(1  + Tilt)  If- 


Expanding  the  (1  + aT  |^)  terms  in  Eq.  (82)  and  collecting 
allows  Eq.  (82)  to  be  rewritten  as 

T w 31  + 111  H *4  - - H *4  + 


terms 


i 4 at  2 "3  at 

2 2 

v.  T 2 t 


T hi. 


M — w if  (St  f 4^7 cos  (qcx)  aSt 


T h T 2 

+ F -1-  COS 


- «c  F ^7  s1n 

- 1c  F 55^-  sin  <V>  W - TF 


_L  ( 

a2T 

ax  \ 

at7 

r 

In. 

3Tc' 

at 

at 

where 


(84) 


W 
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W 
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a2 

(c0 

2 

n n 

9t2 

+ 0 0 
Vo 

92 

- f— 

T h 2 
oo 

It2 

Ip 
\ 0 

3qopo 

(85) 


If  (w0/2tt)  is  the  Fourier  frequency  giving  the  dominant  contri- 
bution to  T,  and  if  the  limit  u)0t  -*•  o (with  t = Tj  or  t2)  is  taken, 

then  Eq.  (83)  may  be  approximated  with 


3 

2 


tiW4 


il 

3t 


Wi 

2Vo 


I 32j  32T 

cos(qcx)  | axat  - J«Si(  + 

I I 

- ic  sin  <V>  ]n-  ir|  • 


(36) 


In  evaluating  the  contributions  from  the  first  and  second  terms  on  the 
right  of  Eq.  (83)  we  have  employed  the  approximation  vQ2  |i-2  ^ 

This  last  approximation  assumes  that  the  temperature  profile  is  of 
the  form 


T(x,t)  = Tj  { q (x  — vQt)  } f (x , t ) 


(87) 


where  qvQ  and  f is  a dimensionless  function  which  varies  slowly 

compared  to  T..  Eq.  (87)  is  dominant  Fourier  term  analysis  in 
disguise,  and*is  conceptually  justified  by  picturing  the  temperature 
profile  as  approximately  in  step  with  the  strain  profile  which  in  turn 
is  propagating  with  a velocity  nearly  v . 


Employing  Eq.  (87)  shows  that  the  ratio  of  the  first  to  second 
term  on  the  right  of  Eq.  (86)  is  given  by  (q/qc)  when  the  signs  and 

cosines  are  taken  as  unity.  Thus  in  the  long  wavelength  limit 
(i.e.  q = 2n/A  goes  to  zero)  we  are  left  with 


W 


H 

4 3t 


- q 


FT  h 
o o 

c 3Vo 


sin  (qcx) 


(88) 
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It  is  possible  to  show  that,  in  the  absence  of  a phase  trans- 
formation (i.e.  with  F=0),  the  strain  satisfies  an  equation  which 
is  one-to-one  with  the  F=0  version  of  Eq.  (83): 


3 

2 


u in  + 
w4  at 


+ 


V T T 
0 12 


(89) 


Thus  in  the  same  (wt)  ■+  0 limit  as  used  in  arriving  at  Eq.  (88)  we 
have 


w4  In  „ n (90) 

at 

in  the  absence  of  transformation  physics.  Comparing  Eq.  (88)  and 
Eq.  (90)  suggests*  that  in  the  presence  of  transformation  effects** 


jjr  = 0 + (second  order  terms).  (91) 

Our  transformation  model,  as  illustrated  in  Fig.  4,  has  the  transition 
temperature  as  a function  only  of  the  strain,  Tc=  Tc(n). 


Thus,  from  Eq. 
3Tc 

u SL  = 

w4  at 


(91), 

in  w 
at  4 


(second  order  terms). 


(92) 


* We  have  not  yet  been  successful  in  deriving  an  explicit  version 
of  Eq.  (91)  with  transformation  effects  included. 

**  If  the  left  hand  side  of  Eq.  (88)  is  considered  to  be  first  order 
(because  of  the  implicit  presence  of  AT  in  aT  ),  then  the  right 

aT 

hand  side  is  second  order  because  of  the  product  F ( the  function 
F contains  the  modal  displacement  uq). 
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If  the  constitutive  law  Tc  = Tc(n)  is  now  chosen  of  the  form 

AT 

= vn  , y = constant,  (93) 

c 

in  keeping  with  experimental  observation14  for  the  room  temperature 
transition  in  TiNi,  then  Eq.  (92)  becomes 


—^C  = 0 + yTc (0>  (second  order  terms),  (94) 

with  T^0)being  the  n = 0 value  of  T . 

Experimentally14  y ^ 10  for  the  TiNi  transition^so  that  one  can 
not,  as  a general  statement,  set  the  right  hand  side  of  Eq.  (94)  to 
zero  compared  to  the  right  hand  side  of  Eq.  (88)  unless  Tj°>  itself 

is  very  small.  We  shall  now  make  that  small  assumption  in  order 

to  facilitate  further  analysis.  This  restriction  to  low  temperature 
transformations  represents  a severe  limitation  on  the  theory  presented 
in  this  report.  Hopefully  the  results  will  at  least  hint  at  the 
behavior  of  room  temperature  transitions. 

The  assumption  of  Eq.  (68)  is  of  a different  class  than  that  of 
small  Tjolin  that  it  is  possible  that  Eq.  (68)  holds  for  some  room 

temperature  transitions. 

-^■y,  see  the  discussion  associated  with  Eq.  (61),  is  a Gruneisen- 
like  parameter.  That  its  magnitude  is  order  ten,  rather  than  of 
order  unity  as  for  non-transforming  media36,  is  probably  very 
significant;  that  same  magnitude  of  approximately  ten  characterizes 
the  Gruneisen  parameter  in  the  region  of  a phase  transition  (see 
Fig.  2).  The  implication  is  that  y is  determined  by  the  physics 
of  the  mode  q . 


32 


Letting  T^°)  be  small  in  Eq.  (94)  gives 


3Tc 

W4^f=  0 


(95) 


so  that  Eq.  (88)  can  be  written 


W4^ 

3t 


- C 


q^  FToho  sin  (q  x)  3ld  , 


3g  p 
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(96) 


Now  set 


T . = T - T 
d c. 


tf2  . 


3Vo  „ 2 
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(97) 


so  that  Eq.  (96)  becomes 


32  I Ha 

It 7 at 


a2 

ax" 


77T  sin  (qcx) 


*4  F 


aT 


_d 

at 


Equation  (98)  is  separable.  Let  ^d  = T (t)x.(x) 

at 

with  Tj  having  the  dimensions  of  , and  Xj  dimensionless 


(98) 


Thus 


a2-, 


4-  4U  = 

^TT  at2 


m2n2 


T " 7" 

'4  '1  ““  '”4  X1 


ax7  ‘ rb  sin  {qcx) 
4 r 


Xj.  (99) 


where  m2n2/c  2 is  the  separation  constant  ( n is  a constant  similar  tow). 

4 0 
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Q2x. 

dx- 


-i  + 


m2q2 


s1n  (V> 

4 F 


X = 0. 

] 


(100) 


The  Laplace  transform,  x^(s),  of  the  function  x(x)  is  defined  by 

xL(s)  =J  x (x)e'SX  dx.  (101) 


Thus  multiplying  both  sides  of  Eq.  (101)  by  e sx  and  integrating 
over  x gives* 


/ o i m2ft2  , / , , i 

(s2  + ) xL(s)  + 22-? 

4 4 r 


xL(s-iqc)  - xL(s+iqc) 


+ sxp)  s «C, 


(102) 


where  the  vertical  slash  with  subscript  zero  indicates  that  the 
functional  to  the  left  is  to  be  evaluated  at  x=0,  and  BC  simply 
denotes  the  boundary  condition  functional. 

If  in  Eq.  (102)  qc  is  set  equal  to  zero  (i.e.  qc=Q)the  result  is 


(s2  + ^ ) x. (s)  = BC  (103) 

which,  upon  taking  the  inverse  transform37  of  X.(s),  is  easily  seen 

to  result  in  propagation  of  d amplitudes  at  constant  phase 

3t 

velocity  (without  growth).  That  same  result  can  be  obtained  directly 
from  Eqs.  (88)  and  (95). 


* The  second  bracketed  term  on  the  left  in  Eq.  (102)  is  obtained 
by  employing  sin(q  x)  = J_exp(iq  x)  - exp(-iq  x)}  prior  to 

c 2i  c c 

integration. 
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Now  consider  zone  edge  qc 
transform  is  given  by37 


(i.e.  qc 


ir/a  in  Fig.  7).  The  inverse 


(104) 


so  that  if  s = ir  we  can  limit  the  range  of  integration  to  |r|  qc 

(from  Eq.  (104)  it  is  seen  that  r plays  the  role  of  q since  r is 
conjugate  to  x,  and  s : ir  transforms  Eq.  (104)  into  a Fourier  trans- 
form if  y is  chosen  to  equal  zero).  Now  let  y = iq  . 


n dnx.  (y) 

:.  x,  (s+iq  ) = x.  (y+s)  = ?_  — r 

L ' x L ' n=o  n!  . n 

dy 


X,  (s-iq)  = X.  (-y+s)  = T.  -j- 


n dnxL(-y) 


n=o 


dy’ 


(105) 


(106) 


Substituting  Eqs.  (105)  and  (106)  into  Eq.  (102)  gives 


(S7  + 
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dy 


xL(y)  - xL(-y) 


sXf (0)  i BC. 


(107) 


If  one  attempts  to  take  the  inverse  transform  of  Eq.  (107)  on  a 
term  by  term  basis,  one  sees37  that  the  contributions  from  n >_  2 in 
the  summation  are  undefined.  Thus,  unless  xL (y ) is  an  even  function 

y,  Eq.  (107)  is  consistent  with  the  possibility  of  avalanching  of  the 

spatial  part  of  8Trf 
u • 

3t 


That  the  qc=0  solution  does  not  predict  avalanching  is  not 
surprising.  qc=0  implies  that  the  transformation  is  associated  with 
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an  infinite  wavelength,  and  consequently  that  the  transformation 
occurs  simultaneously  in  the  entire  solid  sample.  Since  the  physics 
being  discussed  involves  a strain  disturbance  propagating  WITHIN  a 
transformable  medium, there  is  nothing  present  in  the  model  being  used 
which  is  capable  of  exciting  the  qc=0  mode. 

The  above  physics  represents  our  start  on  the  interplay  between 
the  microscopic  and  the  macroscopic  for  strain  pumped  structural 
transformations.  While  the  work  is  very  incomplete  it  is  supportive 
of  the  possibility  of  avalanching  of  the  macroscopic  parameter 
(e.g.  strain  and  temperature). 

Future  work  in  this  area  will  probably  concentrate  on  deriving 
functionally  explicit  versions  of  Eqs.  (91)  and  (94),  and  on  solving 
equations  similar  to  Eq.  (102). 


IV.  THE  a - e TRANSITION  IN  IRON 

Drummond38,  in  1957,  predicted  the  possibility  of  a rarefaction 
(i.e.  expansion)  shock  being  produced  by  the  unloading  process  in  a 
medium  which  has  been  forced  through  a dynamic  (compressive  shock 
induced)  phase  transition.  In  the  absence  of  such  a phase  transition 
(or  other  peculiar  phenomena39*40)  the  unloading  process  whould  pro- 
duce a rarefaction  fan*,41*42  rather  than  a shock. 

The  intersection  of  two  rarefaction  waves  results  in  tension43, 
and  an  adequate  stress-time  history44*45  in  the  tension  regime  can 
cause  spallation  (planar  fracture).  It  has  been  shown  by  Duvall  et  al11 
that  the  relaxation  time  with  which  the  phase  change  induced  mass 
fraction  approaches  an  equilibrium  value  has  an  important  role  in 
determining  the  rise  time  (stress-time  history)  of  the  rarefaction 
shock.  Thus,  since  stress-time  history  is  important  to  spallation, 
it  is  to  be  expected  that  control  of  phase  change  parameters  will 
lead  to  some  control  over  the  satisfying  of  spallation  criteria. 

Such  control  could  be  important  to  the  degree  of  vulnerabi 1 i ty  of 
military  armor. 


* By  fan  we  simply  mean  the  stronger  the  unloading  amplitude,  the 
slower  the  propagation  velocity  associated  with  that  amplitude. 
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Among  the  material  parameters  which  are  important  to  the  Duvall 
et  al  relaxation  time,  percentage  alloy  (e.g.  x in  T i Ni j _x ) and 

impurity  (and  defect)  concentration  are  probably  the  most  significant. 

External  (or  applied)  parameters  probably  also  play  an  important 
role  in  the  phase  change  spallation  problem.  The  a-phase  of  iron  is 
ferromagnetic,  and  one  would  thus  expect  the  presence  of  an  externally 
applied  magnetic  field  to  influence  the  transition.  The  analysis  of 
the  previous  two  sections  involved  the  dominant  Fourier  amplitude 
of  the  applied  strain  disturbance.  The  relative  strength  of  the  strain 
induced  Fourier  amplitudes  is  a function  of  the  rise  time  of  the 
applied  deformation,  and  that  rise  time  can  be  influenced  by  impactor 
structure.  A long  rise  time  can,  for  example,  be  achieved  employing 
an  impactor  of  graded  mass  density  (the  low  mass  density  region  impact- 
ing first). 

At  this  moment  the  details  of  the  above  spall  control  concepts 
are  rather  nebulous.  There  has,  however,  been  some  discussion  in  the 
literature  concerning  the  effect  of  an  applied  magnetic  field  on  the 
a -*  c transition,  and  we  will  now  concentrate  on  the  physics  relevant 
to  that  effect. 

Curran  and  Hornmann46 >47  have  performed  shock  experiments  which 
indicate  that  a magnetic  field  can  reduce  the  dynamic  transition 
pressure  by  approximately  10  kbar.  Barker  and  Hollenbach  subsequently 
performed  a similar  experiment48,  and  found  the  transition  pressure 
unchanged  to  within  approximately  1 kbar.  To  make  the  situation  even 
more  intriguing,  Curran  and  Hornemann46 ,47  have  presented  a derivation 
which  predicts  an  increasing  transition  pressure  with  increasing 
magnetic  field  intensity. 

We  will  now  present  a derivation  which  is  identical  in  its  essentials 
to  that  of  Curran  and  Hornemann  except  that  it  includes  a magnetoelastic49 
contribution. 

Consider  the  Gibb's  free  energy50,  G,  for  a system  composed  initially 
(i.e.  at  atmospheric  pressure)  of  a iron  and  in  which  the  a and  e 
phases  can  exist  in  equilibrium  at  a pressure  P (H)  where  c denotes 

critical  and  H is  the  magnetic  field  intensity. 


G = E - TS  + PV  (108) 
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where  E is  the  total  internal  energy  of  the  system,  and  S the  entropy. 


.•  dG  = dE  - TdS  - SdT  + PdV  + VdP, 


(109) 


but,  assuming  a reversible  process,  from  the  first  law  of  thermodynamics 


TdS  = dE  + PdV  (110) 


so  that  Eg.  (109)  becomes* 


dG  = - SdT  + VdP.  (Ill) 

If  it  is  now  assumed  that  the  temperature  changes  are  insignificant 
as  compared  to  the  pressure  changes,  then 

,P  (H) 

G (PC(H),  T)  - G (0,T)  = / C V(P,T,H)  dP,  (112) 

o 

where  we  include  a V dependence  on  H in  order  to  allow  for  magneto- 
elasticity. The  magnetoelasticity  of  interest  here  arises  due  to  the 
interaction  of  H with  the  ferromagnetic  structure  of  the  a phase. 

The  dependence  of  P on  H only  mirrors  the  fact  that  magnetostriction 

makes  a contribution  to  the  effective  elastic  constants.  Thus  the 
pressure  necessary  to  achieve  a given  critical  volume  is  magnetic 
field  dependent. 


* In  writing  Eqs.  (108)  through  (111)  we  have  implicitly  included 
magnetic  energy,  proportional  to  H , in  E.  We  have  also  assumed 
that  the  magnetic  field  frequencies  of  interest  are  so  small  that 
they  do  not  make  a separate  pressure  contribution  and  influence  V 
through  the  electrical  conductivity  related  penetration  depth. 
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Now  write  Eq.  (112)  separately  for  each  phase. 

G0(pc(H).T)  - G^H)(0,T) 

p (h) 

r c 

= / V (P,T,H)  dP, 
J 1 

(113) 

0 

Ge(Pc(H),T)  - g[H)(0,T) 

fPr(H) 

= J c Ve(P,T,H)  dP, 

(114) 

0 


where  Eq.  (114)  represents  an  extension  of  the  e phase  equation  of 
state  down  to  P=0.  But  a system  of  two  phases  (i.e.  a and  e)  which 
are  in  equilibrium  at  constant  temperature  and  pressure  (i.e.  at 
constant  P (H)  ) has  the  property50 


G (P  (H),T)  = G (P  (H) ,T) . 

a C c C 


(115) 


Thus  upon  subtracting  Eq.  (114)  from  Eq.  (113)  while  using  Eq.  (115)  gives 

(116) 


g*h,«),t)  - g)h'(o,t) 


rVH) 

■J  k-M 


dP. 


-1  0 


-H, 


If  the  magnetic  field  contribution  to  G (0,T)  is  treated  as  a 
perturbation,  then 


G(H)  (0,T)  = G(0)  (0,T)  + g(H) , 


(117) 


where51 


g(H)  = - i HMV  (118) 

In  Eq.  (118)  M denotes  the  magnetic  moment  per  unit  volume,  and  the 
factor  of  h arises  from  assuming  that  M is  linear  in  H.  Experimentally52 
it  is  known  that  the  magnetic  moment  of  the  e phase  is  zero  (i.e.  the  c 
phase  is  not  ferromagnetic)  so  that  g = 0.  Furthermore,  for  H = 0, 
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jGa(0)(0,T)  - G£ (0)(0,T) 


I 


Pr(0) 

Lrv  - v 

1 a e 


H=0 


Thus,  from  Eqs.  (116),  (118),  and  (119) 


-J 


PC(H) 

' ^ H 


dp  - 


rpc(0) 

| c (aV)q  dP. 


dP. 


(119) 


(120) 


Equation  (120)  allows  a solution  for  6P  s Pc ( H ) - Pc ( 0 ) . First 
(aV)h  is  to  be  expanded  to  first  order  in  H. 


3(aV). 

<aV>H  ■ <4V>0  * -V 


with  the  result 


l rPc^H^(AV)H  | 

j HMV  2:  (AV)0  (6P)  + H J — 


dP. 


(121) 


(122) 


If,  in  Eq.  (122), the  partial  derivative  with  respect  to  H is  set 
equal  to  zero,  one  would  essentially  have  the  theoretical  result  of 
Curran  and  Hornemann46’47  with  the  consequence  that  6P  would  be 
predicted  to  be  positive.  Let  6Pt  denote  that  value  of  <5P. 

: <av)0  1 = i 

•SPj  > ~ j »H.  (123) 

<5 P j will  now  be  evaluated  using  the  following  values  for  the  relevant 
parameters  in  Eq.  (123). 
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pm  -j 

Volume  terms53:  v (130  kb,  H=0)  = 6.62  — y- 

a mo  i e 


,e(130  kb,  H-0)  ■ 6.29 


Vo  5 2 (,a  * Vc>-  (lv)0  s <va  • vc>' 


(124) 


(AV) 


- = 5.11  x 10'2 


Field  Intensity46:  H = 2000  Ampere- turns/meter 

Magnetization  (M):  M = KpQH. 


(125) 


In  Mks  units  = 4tt  x 10 


Henrys/meter , and  K is  the  dimensionless 
susceptibility54.  Substituting  the  above  in  Eq.  (123)  gives 


6P  = 

l 


(4.9  x 10  ) K Ne^7~,nS 


(4.9  x 10"7)K  kbar. 


(126) 


For  a iron  K has  an  approximate  range  of  between  120  and  5 x 10J 
(depending  upon  previous  magnetic  history  and  H).  A magnetic  field 
intensity  of  2000  Ampere- turns/m  (equivalent  to  25  Gauss)  is  small  ^ 
compared  to  the  saturation  magnetic  field  for  iron  of  between  2 x 10J 
Gauss  and  4 x 10^  Gauss.  Thus,  while  understanding  that  previous 
magnetic  history  is  most  important  for  small  values  of  H,  it  seems 
most  reasonable  to  use  K = 2 x 10^  in  Eq.  (126),  with  the  result 


6Pj  -v  10" 


kbar. 


(127) 


The  e phase  is  not  ferromagnetic,  and  its  magnetostriction  can 
thus  be  neglected.  Additionally  we  will  assume  that  the  a phase 
magnetostriction  is  independent  of  stress,  with  the  result  that  Eq . ( 1 22 ) 
can  be  approximated  by 


(SP) 


HP 


3V 
c 

3H 


(128) 
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where  for  approximation 
The  problem  now  reduces 


purposes  we  have  replaced  P (H)  by  P =130  kbar. 
3V  C C 

to  evaluating  _a  from  magnetostriction  data. 
8H  o 
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Figure  8.  Typical  experimental 
arrangement. 


Figure  8 illustrates  the  typical 
arrangement  of  the  type  of  experiment 
being  discussed  here.  The  field  H 
can  be  induced  by  a permanent  magnet 
or  by  having  the  target  be  an  integral 
part  of  an  electromagnet.  If  it  is 
imagined  that  the  magnetostricti ve 
effect  is  associated  with  a stress 
in  the  direction  of  M in  Fig.  8,  then 


(129) 


where  v is  Poisson's  ratio20,  and 
(6t/t)  is  the  strain  in  the  M direction.  For  brevity  let  v = 1/3. 
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10  L 


Figure  9 crudely  illustrates 
magnetostriction  data  for  iron 
(from  Mattiat49).  For  the  2000 
Ampere- turn/m  external  field  value 
of  Curran  and  Hornemann  we  thus 
have,  upon  using  Eq.  (129), 


\ 

\ 

5-f- 


-il 


Tto" 


H 

TT 

o 


3 V 

a 

w 


=W6- 


(130) 


^2 


H x 10 

(A-turns/m)  Letting  SPr.  be  the  contribution 

~2to"  t0  ^rom  the<i[Tia9net0Str''ct'' ve 


effect 


\ 

/ 
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HP 

3V  1 

6P2,  -( 

0 

c 
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^ol 

Vo 

3H 

(131) 


Fiqure  9.  A crude  approximation 
to  the  magnetostriction 
data  for  iron. 
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T 


«p2  * - SA  ;■  }S-‘  PC  ' - 5-5  * I0'3  kb-  <»2> 


The  most  interesting  aspect  of  of  Eq.  (132)  is  its  sign;  the 

sign  of  6P^  indicates  a decrease  in  transition  pressure  in  the  presence 

of  a magnetic  field,  in  agreement  with  the  experiments  of  Curran  and 
Hornemann.  Of  course  Curran  and  Hornemann  observed  a decrease  in 
transition  pressure  of  approximately  10  kbar,  while  our  6P£  value  is 

three  orders  of  magnitude  smaller  than  the  mentioned  experimental 
observations. 

The  most  drastic  assumption  made  in  arriving  at  Eq.  (132)  is 
probably  that  of  stress  independent  mangetostriction.  The  applied 
magnetic  field  as  well  as  the  shock  wave  each  have  an  orientation 
effect  upon  the  magnetic  domain  structure  in  an  iron  sample.  While 
it  is  difficult  to  imagine  how  such  an  orientation  interaction  can 
effect  6 by  three  orders  of  magnitude,  that  possibility  should 

not  be  immediately  ruled  out. 

Barker  and  Hoi lenbach' s48  experiments  were  carried  out  at  with 
an  external  field  of  H = 2000  Gauss,  and  observed  SP  = 0 to  within 
1 kbar.  But  2000  Gauss  is  equal  to  1.6  x 10$  Ampere-turns  per  meter, 
and  thus  represents  a magnetic  field  approximately  two  orders  of 
magnitude  larger  than  that  employed  by  Curra"  and  Hornemann.  It  is 
most  interesting  to  note  that,  at  least  as  far  as  the  magnetostriction 
effect  is  concerned  (see  Fig.  9),  the  different  magnetic  field 
magnitudes  give  potentially  radically  different  experiments.  It  thus 
seems  that  further  experiment  and  theory  must  be  carried  out  before 
this  issue  can  be  resolved.  There  is  an  applied  motivation  for  such 
further  work;  a 6P  >10  kbar  could  play  a role  in  the  further  engi- 
neering development  of  tank  armor;  for  a significant  6P  one  could 
envision  a tank  commander  flipping  a switch  in  order  to  magnetize 
his  vehicle  just  prior  to  an  expected  fragment  impact. 


V.  CONCLUDING  DISCUSSION 

This  report  has  been  concerned  with  the  interplay  between  phase 
change  effects  and  shock  wave  physics.  The  report  represents  a 
summary  of  some  beginning  ideas  on  the  subject. 

Two  motivated  problems  have  been  of  main  concern.  The  study  of 
shock  induced  structural  phase  transitions  in  metallic  alloys  was 
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motivated  by  a desire  to  obtain  a more  efficient  conversion  medium 
for  use  in  an  optically  activated  detonator.  The  consideration  of 
the  effect  of  an  applied  external  magnetic  field  on  the  a -*•  e phase 
transition  in  iron  was  motivated  by  a desire  to  investigate  the  role 
of  some  esoteric  solid  state  concepts  on  the  vulnerability  of  armor. 

for  the  optically  activated  detonator  problem  our  specific  interest 
was  in  the  possibility  that  the  shock  induced  transition  might  result 
in  avalanching  of  the  amplitude  of  the  propagating  shock  wave.  Both 
macroscopic  and  microscopic  considerations  are  shown  to  be  consistent 
with  the  possible  existence  of  such  avalanching.  In  order  to  reach 
more  definite  conclusions  further  work  is  necessary.  That  further 
work  includes  more  rigorous  modeling  of  the  macro  and  micro  aspects 
of  transition  physics  as  well  as  the  seeking  of  solutions  of  equations 
which  are  in  the  form  of  Eq.  (102).  The  more  rigorous  modeling  includes 
replacing  the  assumption  of  Eq.  (68)  with  better  physics,  and  attempt- 
ing to  explicitly  include  shear  (wave)  aspects  of  the  transformation. 

The  displacements  associated  with  the  room  temperature  transition  in 
NiTi  correspond  to  shear  modes,  and  the  one  dimensional  approach 
presented  in  this  report  has  ignored  that  very  important  fact. 

Although  the  report  predicts  a very  small  magnetoelastic  contri- 
bution (of  the  order  of  10“2  kbar)  to  the  a ->■  e transition  pressure 
in  iron,  that  prediction  is  very  preliminary.  Considerably  more 
work  remains  to  be  done  in  this  problem  area,  for  both  magnetic 
field  influences,  and  with  respect  to  the  possibility  of  tweaking 
other  materials  parameters  which  can  be  important  to  transformation 
physics. 


44 


REFERENCES 


1.  J.  A.  Krumhansl  and  J.  R.  Schrieffer,  Phys.  Rev.  B 1_1,  3535  (1975). 

2.  T.  Schneider  and  E.  Stoll,  Phys.  Rev.  B 1J3 , 1216  (1976). 

3.  See,  for  example,  G.  B.  Whitham,  Unear  and  Nonlinear  Waves 
(John  Wiley  and  Sons,  New  York,  1974). 

4.  L.  C.  Yang  et  al.  National  Defense,  page  344  (January  - February  1974). 

5.  L.  C.  Yang,  J.  Appl . Phys.  45,  2601  (1974). 

6.  F.  E.  Walker  and  R.  J.  Wasley,  Explosivstoffe  17,  9 (1969). 

7.  N.  G.  Pace  and  G.  A.  Saunders,  Solid  State  Comm.  9,  331  (1971). 

8.  For  a thorough  review  of  the  properties  of  TiNi  alloys  see 

C.  M.  Jackson  et  al,  NASA  Technical  Report  NASA-SP  5110  (1972) 

Available  from  Superintendent  of  Documents,  U.S.  Government 
Printing  Office,  Washington,  DC  20402. 

9.  This  is  not  the  first  time  that  the  US  Army  has  shown  interest 
in  NiTi  alloys  for  fuze  train  applications.  To  the  best  of  our 
knowledge,  however,  all  of  the  other  studies  have  been  concerned 
with  the  "shape-memory"  effect.  See,  for  example,  C.  Friedman, 
Frankford  Arsenal  Memorandum  Report  M68-38-2,  November,  (1968). 
Frankford  Arsenal,  Philadelphia,  PA  19137. 

10.  F.  E.  Wang  et  al , J.  Appl.  Phys.  36,  3232  (1965). 

11.  G.  E.  Duvall  et  al,  "Equation  of  State  of  Solids  (Final  Report), 
Ballistic  Research  Laboratories  Contract  Report  No.  67,  May  1972. 
Prepared  by  Washington  State  Univ.,  Pullman,  Washington,  as 
report  WSU-SDL  71-01. 

12.  W.  J.  Murri  et  al,  article  in  Advances  in  High  Pressure  Research, 
vol . 4,  R.  H.  Wentorf,  Jr.,  editor  TAcademic  Press,  New  York,  1974). 

13.  M.  W.  Zemansky,  Heat  and  Thermodynamics,  third  edition  (McGraw- 
Hill,  New  York,  1951) . 

14.  See,  for  example,  figure  78  of  reference  8. 

15.  A.  D.  Smirnov,  Tables  of  Airy  Functions  and  Special  Confluent 
Hypergeometric  Functions  (Perqamon,  New  York,  1960). 


45 


16.  The  results  for  Cj  = c J2  are  to  appear  in  a paper  by  P.  Harris 

in  the  Transactions  of  the  22nd  Conference  of  Army  Mathematicians, 
Watervliet  Arsenal,  Hay  1976  (Army  Research  Office,  Durham,  North 
Carol ina). 

17.  A.  Jeffrey,  Int.  J.  Non-Linear  Mech.  6,  669  (1971). 

18.  R.  F.  Bunshah  and  R.  F.  Mehl , Trans.  AIMt  - J.  Metals,  page  1251, 
September  1953. 

19.  W.  J.  Moore,  Seven  Solid  States  (Benjamin,  New  York,  1967). 

20.  H.  Kolsky,  Stress  Waves  in  Solids  (Dover,  New  York,  1963). 

21.  J.  M.  Ziman,  Electrons  and  Phonons  (Clarendon  Press,  Oxford, 
England,  1962T! 

22.  A.  A.  Maradudin,  E.  W.  Montroll,  and  G.  H.  Weiss,  Theory  of 
Lattice  Dynamics  in  the  Harmonic  Approximation  (Academic  Press, 

New  York,  1963).  This  monograph  is  supplement  #3  of  the  Solid 
State  Physics  series. 

23.  W.  Cochran,  Reports  on  Progress  in  Physics  26,  1 (1963). 

24.  J.  deLaunay  in  Solid  State  Physics  (F.  Seitz  and  D.  Turnbull,  eds). 
Vol . 2.  (Academic  Press,  New  York,  1956). 

v *i 

25.  R.  Blinc  and  B.  Zeks,  Soft  Modes  in  Ferroelectrics  and  Anti- 
ferroel ectrics  (American  Elsevier,  New  York,  T974). 

26.  J.  D.  Axe  and  G.  Shirane,  Physics  Today,  Sept.  1973,  page  32. 

27.  T.  Schneider,  G.  Srinivasan,  and  C.  P.  Enz,  Phys.  Rev.  A 
1528  (1972). 

28.  For  example  see  the  illustration  on  page  34  of  reference  26. 

29.  L.  D.  Landau  and  E.  M.  Lifshitz,  Statistical  Physics  (Addison- 
Wesley,  Reading,  1958).  Section  134. 

30.  K.  R.  Atkins,  Liquid  Helium  (Cambridge  University  Press, 

Cambridge,  1959) . 

31.  T.  H.  K.  Barron,  Phil,  Mag.  46,  720  (1955). 

32.  K.  Brugger,  Phys.  Rev.  137,  A1826  (1965). 


46 


33.  T.  0.  Woodruff  and  H.  Ehrenreich,  Phys.  Rev.  1^3,  1553  (1961). 

34.  H.  H.  Barrett,  Phys.  Rev.  178,  743  (1969). 

35.  P.  Harris,  Technical  Report  4902,  Picatinny  Arsenal,  Dover,  NJ, 
October  1975. 

36.  P.  Harris  and  L.  Avrami,  Technical  Report  4423,  Picatinny  Arsenal, 
Dover,  NJ  07801,  Sept.  1972. 

37.  G.  E.  Roberts  and  H.  Kaufman,  Table  of  Laplace  Transforms 
(W.  B.  Saunders  Co.,  Philadelphia,  1966). 

38.  W.  E.  Drummond,  J.  Appl . Phys.  28,  993  (1957). 

39.  P.  Harris,  J.  Acoust.  Soc.  Amer.  4J1,  226  (1966). 

40.  L.  M.  Barker  and  R.  E.  Hollenbach,  J.  Appl.  Phys.  4U  4208  (1970). 

41.  R.  Courant  and  K.  0.  Friedrichs,  Supersonic  Flow  and  Shock  Waves 
(Interscience,  New  York,  1948).  Section  45. 

42.  Ya.  B.  Zeldovich  and  Yu.  P.  Raizer,  Physics  of  Shock  Waves  and 
High  Temperature  Hydrodynamic  Phenomena  (Academic  Press,  New  York, 
1966).  Section  17  of  Chapter  I,  and  Section  20  of  Chapter  XI. 

43.  B.  M.  Butcher  et  al,  AIAA  Journal  2,  977  (1964). 

44.  F.  R.  Tuler  and  B.  M.  Butcher,  Int.  J.  of  Frac  Mech.  4,  431  (1968). 

45.  J.  C.  Peck,  H.  M.  Berkowitz,  and  L.  J.  Cohen,  Int.  J.  Frac.  Mech. 

5,  297  (1969). 

46.  D.  R.  Curran  and  U.  Hornemann,  "Der  Einfluss  eines  Magnetfeldes 
auf  die  Hochdruckphasendiagramme  von  Armco-Eisen  und  normiertem 
Baustahl,"  Bericht  Nr.  6/70,  Ernst-Mach-Institut,  Freiburg/Br, 

West  Germany  (1970). 

47.  D.  R.  Curran,  in  Shock  Waves  and  the  Mechanical  Properties  of 
Sol  ids,  edited  by  J.  J.  Burt-e  and  V.  Weiss  (Syracuse  University 
Press,  Syracuse,  1971). 

48.  L.  M.  Barker  and  R.  E.  Hollenbach,  J.  Appl.  Phys.  45,  4872  (1974). 

49.  0.  E.  Mattiat,  Ultrasonic  Transducer  Materials  (Plenum  Press, 

New  York,  1971). 


47 


50.  F.  Reif,  Statistical  and  Thermal  Physics  (McGraw-Hill,  New  York,  1965). 

51.  See  Section  11.1  of  reference  50. 

52.  R.  N.  Keeler  and  A.  C.  Mitchell,  Solid  State  Comni.  7,  271  (1969). 

53.  See  figure  3.8  of  reference  11. 

54.  W.  K.  H.  Panofsky  and  M.  Phillips,  Classical  Electricity  and 
Magnet i sm  (Addison-Wesley,  Reading,  1955).  Chapter  8, 


48 


DISTRIBUTION  LIST 


Copy  No. 

Department  of  the  Army 

Office,  Chief  of  Research  and  Development 

ATTN:  Dr.  J.  I.  Bryant  1 

Washington,  DC  20438 

Vice  President 

Sandia  Laboratories,  Livermore 

ATTN:  G.  R.  Otey,  3157  2 

R.  S.  Jacobson  3 

Livermore,  CA  94550 

Di rector 

Lawrence  Livermore  Laboratory 

ATTN:  Dr.  Dick  Weingart  4 

Dr.  Frank  E.  Walker  5 

Livermore,  CA  94550 

Kaman  Sciences 
Garden  of  the  Gods  Road 

ATTN:  C.  W.  Gullikson  6 

D.  Williams  7 

Colorado  Springs,  CO  80907 

Effects  Technology  Corporation 

ATTN:  F.  Tuler  8 

5383  Hollister  Avenue 

P.0.  Box  30400 

Santa  Barbara,  CA  93105 

Stanford  Research  Institute 
Poulter  Laboratories 

ATTN:  Dr.  William  J.  Murri  9 

Dr.  D.  R.  Curran  10 

Dr.  R.  K.  Linde  11 

Menlo  Park,  CA  94025 

Commander 

Harry  Diamond  Laboratories 

ATTN:  Mr.  Philip  Brody  12 

Dr.  R.  B.  Oswald,  Jr.  13 

Washington,  DC  20438 


49 


14 


Commander 

USAECOM,  CS&TA  Laboratory 

ATTN:  AMSEX-CT-L,  Dr.  R.  G.  Buser 

Fort  Monmouth,  NJ  07703 


Sandia  Corporation 

ATTN:  Dr.  Walter  Hermann  lb 

Dr.  Robert  Graham  16 

Dr.  J.  C.  King  (1900)  17 

Dr.  D.  H.  Anderson  (1910)  18 

Dr.  J.  Gover  (1935)  19 

Dr.  William  Benedick  20 

Dr.  R.  E.  Hollenbach  21 

Dr.  L.  D.  Bertholf  22 

P.0.  Box  5800 

Albuquerque,  NM  87116 

Terra  Tek,  Inc. 

ATTN:  Dr.  L.  M.  Barker  23 

420  Wakara  Way 

Salt  Lake  City,  UT  84108 

Washington  State  University 

ATTN:  Dr.  George  Duvall  24 

Dr.  George  Swan  25 

Pullman,  WA  99163 

Dr.  D.  V.  Keller,  President 

KTech  Corp.  26 

911  Pennsylvania  Ave.,  N.E. 

Albuquerque,  NM  87110 


Commander 

U.S.  Naval  Ordnance  Laboratory 
Explosion  Dynamics  Division 

ATTN:  Dr.  D.  John  Pastine  27 

Dr.  S.  J.  Jacobs  28 

Dr.  J.  Forbes  29 

Dr.  James  Goff  30 

White  Oak,  Silver  Spring,  MD  20910 

Comma nder 

US  Army  Research  Office 

ATTN:  Dr.  E.  A.  Saibel  31 

Dr.  J.  Chandra  32 

P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 


50 


Commander 

US  Army  Research  & Standardization  Group  (Europe) 

ATTN:  Dr.  A.  K.  Nedoluha  33 

COL  J.  M.  Kennedy,  Jr.  34 

P.0.  Box  65 
FPO  09510 

National  Bureau  of  Standards 

ATTN:  Dr.  Donald  Tsai  35 

Dr.  Henry  Prask  36 

Gaithersburg,  MU  20760 

California  Institute  of  Technology 

ATTN:  Dr.  Thomas  J.  Ahrens  37 

Pasadena,  CA  91109 

Commander 

Ballistic  Research  Laboratories 

ATTN  Dr.  Stan  M.  Taylor  38 

Dr.  George  Adams  39 

Dr.  Robert  F.  Eichelberger  40 

Mr.  George  E.  Hauver  41 

Dr.  Denis  F.  Strenzwilk  42 

Dr.  I.  May  43 

Aberdeen  Proving  Ground,  MD  21005 

University  of  Illinois 

Department  of  Chemistry  and  Chemical  Engineering 

ATTN:  Dr.  H.  G.  Drickamer  44 

Urbana,  IL  60436 

Commander 
Watervliet  Arsenal 

ATTN:  Dr.  Y.  K.  Huang  . 45 

Dr.  T.  E.  Davidson  46 

Watervliet,  NY  21289 

Physics  International  Company 

ATTN:  Dr.  John  Huntington  47 

Dr.  James  Shea  48 

2700  Merced  Street 
San  Leandro,  CA  94577 

University  of  Delaware 
Department  of  Physics 

ATTN  Prof.  Ferd  E.  Williams  49 

Prof.  W.  B.  Daniels  50 

Newark,  DE  19711 


51 


Director 

Defense  Documentation  Center 
ATTN:  DDC-TCA 

Cameron  Station,  Building  5 
Alexandria,  V A 22314 

Union  Carbide  Corporation 
Tarrytown  Technical  Center 
ATTN:  Dr.  John  B.  Lightstone 

Dr.  Jaak  Van  Den  Sype 
Tarrytown,  NY  10591 

Commanding  General 
US  Army  Materiel  Command 
ATTN:  AMCRD-X,  Dr.  Halley 

Washington,  DC  20315 

McDonnel  Douglas  Astronautics 
ATTN:  Dr.  John  Watcher 

5301  Bolsa  Ave. 

Huntington  Beach,  CA  92647 

Systems,  Science,  and  Software 
ATTN:  Dr.  H.  E.  Read 

P.0.  Box  1620 
La  Jolla,  CA  92037 

Director 

Defense  Nuclear  Agency 

ATTN:  SPAS,  Mr.  J.  F.  Moulton,  Jr. 

Washington,  DC  20305 

Army  Materials  & Mechanics  Center 
ATTN:  Mr.  John  F.  Dignam 

Mr.  John  Mescal  1 
Bldg.  131 
Arsenal  Street 
Watertown,  MA  02172 

Lockheed  Palo  Alto  Research  Labs 
ATTN  Dr.  J.  F.  Riley 
3251  Hanover  St. 

Palo  Alto,  CA  94304 

Director 

Los  Alamos  Scientific  Laboratory 
ATTN:  Dr.  J.  M.  Walsh 

Los  Alamos,  NM  87544 


51-61,  142 


62 

63 


64 

65 


66 


67 

68 
69 


70 


71 


52 


T 


North  Carolina  State  University 

ATTN:  Prof.  Y.  B.  Horie 

Dept,  of  Engineering  Science  & Mechanics 

Raleigh,  NC  27607 

University  of  Tennessee 
ATTN:  Prof.  M.  A.  Breazeale 

Dept,  of  Physics  and  Astronomy 
Knoxville,  TN  37916 

Director 

National  Bureau  of  Standards 

ATTN:  Dr.  V.  Arp,  Cryogenics  Division 

Boulder,  CO  80302 

Mellon  Institute 

ATTN:  Dr.  Bernard  D.  Coleman 

Pittsburgh,  PA  15213 

Princeton  University 
ATTN:  Prof.  A.  C.  Eringen 

Prof.  Peter  Mark 
Princeton,  NJ  08 540 

Carnegie  Institute  of  Technology 
ATTN:  Prof.  Morton  E.  Gurtin 

Dept,  of  Mathematics 
Pittsburgh,  PA  15213 

Brown  University 
ATTN:  Prof.  Robert  T.  Beyer 

Dept,  of  Physics 
Providence,  RI  02912 

Courant  Institute  of  Mathematical  Sciences 

ATTN:  Library 

New  York  University 

New  York,  NY  10453 

Queens  College  of  the  City  University  of  New  York 
ATTN:  Prof.  Arthur  Paskin 

Dept,  of  Physics 
Flushing,  NY  11300 


72 


73 


74 


7b 


76 

77 


73 


79 


80 


81 


53 


Commander 
Picatinny  Arsenal 

ATTN:  SARPA-CO,  Mr.  H.  W.  Painter  82 

Chief,  Technical  Information  Br  83-92 

SARPA-FR,  Dr.  E.  Sharkoff  93 

SARPA-FR-E,  Dr.  R.  Walker  94 

SARPA-FR-E-S,  Dr.  Harry  Fair  95 

SARPA-FR-E-S,  Mr.  Louis  Avrami  96 

SARPA-FR-E-S,  Dr.  T.  Gora  97 

SARPA-FR-E-S,  Dr.  F.  Owens  98 

SARPA-FR-E-S,  Dr.  Igbal  99 

SARPA-ND,  Mr.  W.  Benson  100 

SARPA-AD,  Mr.  A.  Moss  101 

SARPA-ND-C,  Mr.  G.  I.  Jackman  102 

SARPA-ND- C- SC,  Mr.  W.  Reiner  103 

SARPA-ND-C,  Dr.  Paul  Harris  104-134 

SARPA-ND-C,  Mr.  A.  Garcia  135 

SARPA-FR-S,  Mr.  J.  Pearson  136 

SARPA-AD-EP4,  Mr.  D.  Seeger  137 

SARPA-FR-S-P,  Mr.  W.  Doremus  138 

SARPA-FR-S-P,  Dr.  G.  Vezzolli  139 

SARPA-ND-C-TS,  Mr.  H.  Opat  140 

SARPA-MI-M,  Mr.  B.  Barnett  141 

Dover,  NJ  07801 


54 


